Valence Bond Glass Phase in the Diluted Kagome Antiferromagnets 
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We present a theory for site dilution in the Valence Bond Crystal Phase of the Kagome Lattice 
Heisenberg Model. The presence of an empty site leads to strong singlet bonds across the impurity. 
It also creates a free spin, which delocalizes inside the unit cell. Finite concentration of quenched 
impurities leads to a Valence Bond Glass phase. This phase has short-range Valence Bond order, no 
spin-gap, large spin-susceptibilities, linear specific heat due to two-level systems, as well as singlet 
and triplet excitations that decompose into kink-anikink pairs delocalized over a few lattice con- 
stants. It provides a framework for a comprehensive understanding of thermodynamic, neutron, and 
Raman measurements in the Herbertsmithite material ZnCu3(OH)6Cl2, including recently reported 
H/T and lj/T scahng. 

PACS numbers: 



The Kagome Lattice Heisenberg Model (KLHM) is one 
of the most studied realistic quantum spin-models where 
Resonating Valence Bond (RVB) physics, first proposed 
by Anderson [l| more than twenty five years ago, may 
be reaUzed. In this model, extensive degeneracy at the 
classical level and in the space of Valence Bond config- 
urations leads to the possibility of quantum spin liquids 
as well as many competing or coexisting orders. j2ri 16]. 

Experimental studies of the Herbertsmithite materials 
ZnCu.-H (011)^0 12, containing structurally perfect Kagome 
planes, 17h23| have raised hope of realizing RVB physics 
in a real material. Ideally, these materials have a py- 
rochlore structure, where spin-half copper atoms form 
Kagome layers, which are separated by non-magnetic 
zinc containing triangular layers. This ideal case is now 
known to be not true as substitution of some fraction of 
zinc and copper sites leads to extra isolated spins in the 
zinc planes and site-dilution in the copper planes. 

Here, we develop a theory for site dilution in KLHM 




in context of Dimer series expansions. [10|, |2J] We find 
that dilution leads to a Valence Bond Glass phase, which 
has short-range Valence Bond order, no spin-gap, large 
spin-susceptibility, linear specific heat due to two-level 
systems as well as singlet and triplet excitations that 
decompose into kink-antikink pairs with energies spread 
upto approximately 2.5J. Unlike the VBC phase, which 
is stabilized by higher order quantum fluctuations, with 
energy scales less than 0.01 J per site (energy scales that 
differentiate different VBC phases are only 0.001 J per 
site), the VBG phase is stabilized by second order per- 
turbation theory with an energy scale of order 0.1 J per 
site. This VBG phase provides a novel framework for 
a comprehensive understanding of the Herbertsmithite 
materials. 

A single impurity can be accomodated in the Valence 
Bond Crystal phase in several inequivalent positions of 
the 36-site unit cell with nearly equal energy. The ad- 
dition of such an impurity leads to strong singlet bonds 
across from it, in agreement with earlier studies of Dom- 



FIG. 1: Proposed Valence Bond Phase of the Kagome Lattice 
Heisenberg Model consists of a Honeycomb Lattice of resonat- 
ing hexagons (H), where six hexagons surround a pinwheel 
(P). The pinwheels are isolated from empty triangles leading 
to substantially reduced quantum fluctuations and doubly de- 
generate dimer configurations. 



mange et al.[25| However, contrary to the work of Dom- 
mange et al, we find that the impurity also creates a free 
spin, which delocalizes inside the unit cell. The Herbert- 
smithite materials have a large impurity concentration of 
order six percent, which would imply on average 2 ran- 
domly placed impurities per unit cell. Combined with the 
singlet bonds across from the impurities, it also means 30 
percent of spins are inactive. At this high impurity con- 
centration no semblence of long range VBC order would 
remain. Instead, one obtains a randomly pinned Valence 
Bond Glass phase with only short range Valence Bond 
order. We will show that the VBG description leads to a 
consistent picture for thermodynamics, neutron and Ra- 
man scattering experiments in the Herbertsmithite ma- 
terials. 

The dimer expansions provide a hirearchy of energy 



rde r q uantum fluctuations in the 
1).[10| A dimer configuration has 



scales for different order 
KLHM ( we set J 
energy —0.375. After including second order quantum 
fluctuations, the energy of each dimer configuration is 
lowered to —0.421875. We emphasize that at this level 
all dimer configurations lead to the same energy. On the 
other hand the exact ground state energy of KLHM is 
known from a variety of numerical studies to be —0.433. 
Thus the entire set of higher order quantum fluctua- 
tions beyond second order change the energy of the sys- 
tem by only another —0.011. The first selection in the 
space of dimers comes from the occurance of resonating 
hexagons. [?] In 3rd order, having a maximally allowed 
number of such hexagons lowers the energy of the state 
to —0.42578125. These hexagons also contribute a frac- 
tion of the fourth order perturbation theory, that brings 
the energy to approximately —0.431. The favored VBC 
configuration, which is selected in fourth order perturba- 
tion theory is shown in Fig. 1. We should note that the 
same VBC configuration has been obtained variationally 
in a completely unbiased tensor-network formalism by 
Evenbly and Vidal.^ 

The addition of holes lifts the degeneracy of the dimer 
configurations already in second order perturbation the- 
ory. For a finite concentration of holes, not having a 
strong bond in a triangle with a missing spin raises the 
energy of the dimer configuration by about 0.05 if one 
bond opposite the missing spin is absent and 0.1 if both 
bonds opposite the missing spin are absent. Given this 
relatively large energy scale, it is fair to assume that these 
impurities freeze all quantum fluctuations except those 
involving (i) a pair of neighboring dimers (second or- 
der perturbation), (ii) resonating hexagons and (iii) pin- 
wheels and analogous structures. Because the pinwheels 
have a doubly degenerate dimer configuration, they play 
a vital role in the VBG phase at low impurity concentra- 
tion. 

We first ask the question, how would the crystal re- 
spond to a single impurity? Given the many inequiva- 
lent sites, the crystal is free to adjust itself to place the 
impurity at any site. Second order perturbation theory 
tells us that the impurity will prefer to go to one of the 
bonds that do not touch any empty triangles. These 
are shown as the dark bonds in Fig. 1. They are the 
same bonds that lead to heavy (non-mobile) triplets in 
the VBC phase. [24*1 To a high numerical approximation, 
there are three inequivalent sites where the missing spin 
can go. These are (i) the inside hexagon of the pinwheel, 
(ii) The outside vertices of the pinwheel and (iii) on one 
of the dark bonds that are not part of the pinwheel. The 
structures are all bounded from outside by dimerized tri- 
angles, which are in their ground state. Hence, it suffices 
to diagonalize the pinwheels only, plus the extra triangle 
in case (iii). 

The ground state bond energies and spin configura- 
tions for the three cases are shown in Fig. 2. Fig. 2(a) 






FIG. 2: Dimer configurations around the impurity, when the 
impurity is placed (i) in the inner hexagon, (ii) on the outer 
vertices, and (iii) in the bond neighboring the pinwheel. Fig 
(a) through (c) show the dimerization pattern, with the thick- 
ness of the lines showing the strength of the singlet bonds. 
Note that the other side of the impurity in (b) and (c) is a 
non-fluctuating singlet bond. Fig (d) through (f) show the 
distribution of free spins. The filled circles show the positive 
spin expectation values, whereas empty circles show negative 
values. The size of the circle represents the size of the spin. 



through (c) show the pattern of dimerization around the 
impurity. Fig. 2(d) through (f) show how the free spin 
is delocalized over the pinwheel. In all cases, the energy 
cost of removing a spin is not 0.75 as for an isolated dimer 
but rather approximately 0.25. The system chooses to 
put a strong singlet bond across the impurity, which gets 
back the lost 0.75 of energy. Instead an antikink devel- 
ops in the pinwheel, which is known to have a minimum 
energy of approximately 0.25. The actual energy cost in 
the three examples is found by exact diagonalization to 
be approximately 0.2764, 0.2632 and 0.2468 in cases (i), 
(ii) and (iii) respectively. The singlet bond across the 
impurity has strength —0.725, —0.716 and —0.731 re- 
spectively. The bond patterns have a clear resemblence 
to the study of Dommange et al. 25 1 

In contrast to the dark bonds, if the impurity is placed 
on the grey bonds, the resulting free spin can delocal- 
izc through the network of grey bonds. It cannot delo- 
calize too far because of the confining dimerizing field, 
but the dimerization between the weak and strong bonds 
connecting the empty triangles is relatively weak. This 
implies that the confining potential is weak and the spin 
can delocalize over several unit cells. 
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With a finite concentration of randomly placed impuri- 
ties, there will be lots of free spins created. When two of 
these spins meet they would bind into a singlet. [31 The 
mobile spins can also wander into the pinwheel regions 
and form singlets with spins there. The strength of the 
pairing will depend on the ability of the spins to be in 
overlapping regions and thus would decrease rapidly with 
separation of impurities. In a thermodynamic system of 
randomly placed impurities, most of the spins would be 
paired into singlets at T near zero. Only very rare re- 
gions with isolated impurities will have a free spin. This 
implies a zero spin-gap for the system. At infinitesimally 
small impurity concentration there would be a 1/T sus- 
ceptibility with a strength proportional to the impurity 
concentration. But, at any significant impurity concen- 
tration this strength would be reduced. We show in ap- 
pendix that the assumption that these spins go into a 
random singlet phase |26| with a power-law distribution 
of exchange constants, allows us to reproduce the H/T 
and uj/T scaling reported recently. [27 

The glassy system has many potential two-level system 
type local excitations created by the impurities. For very 
dilute impurities, the pin- wheels would serve as two level 
systems. The random environment would produce a weak 
random exchange field on the pinwheels giving rise to a 
small splitting, between its two ground states. In a higher 
concentration range, as relevant for the Herbertsmithite 
materials, very few intact pinwheels are likely to remain. 
Instead, at larger impurity concentration, the two-level 
systems would be a network of corner sharing triangles, 
where the Valence Bonds can switch between alternate 
dimer configurations, with very little energy cost. These 
would resemble Delta chains or Husimi trees, 15, 35 



a pinwheel being just one special case. If we assume one 
two level system per 36-sites and further assume that 
their splitting can range upto 10 percent of J (approx- 
imately 20 K) that would lead to a density of states of 
1/{12{)K) per copper atom. This will give rise to a lin- 
ear term in specific heat of approximately 20mJ/molK^ 
per copper atom. This is in the right ballpark of the 
experimental observations on the Herbertsmithite mate- 
rials when the spin degrees of freedom associated with 
free spins in the zinc planes, are supressed by a magnetic 
field. [28] 

The very low energy magnetic response will depend 
on the DM interactions, 2^31| which are not considered 
here. For example, the NMR spectra has already been 
addressed by Rousochatzakis et al in context of frozen 
dimers and it requires non-zero DM interactions in addi- 
tion to impurities. [3^, |33| The VBG phase also provides 
a natural explanation for the observed Neutron spectra 
in the materials. One of the key observations of recent 
powder diffraction neutron measurement of the materi- 
als by deVries et al 3J| is that the spectra is spread over 
a large frequency range and its angle averaged behavior 
is very close to that of isolated dimers. This is exactly 





FIG. 3: Kink-antikink or 2-spinon states for the pinwheel. A 
triplet excitation is created by breaking a singlet bond in the 
ground state. The Kink spinon lies on the inner hexagon, 
while the more mobile antikink spinon can move on the out- 
side vertices of the pinwheel through states (a) through (f) 
shown in the figure. 



what one expects for pinwheels or Delta chains, as we 
show below. 

Because the ground state of a pinwheel is fully 
dimerized, the equal-time correlation function is strictly 
nearest-neighbor only, that is, it is that of isolated 
dimers. However, the triplet excitations in these systems 
decompose into a kink-antikink pair as shown in Fig. 3. 
Hence the dynamic correlations arc not strictly nearest- 
neighbor only and unlike a single dimer they extend over 
a wide energy range. It is known 35, [36| that the kinks 
are gapless and anti-kink energy can be approximated by 



e{k) — 5/4 — cos k, 

so that the kink antikink pair energy ranges approxi- 
mately from about 1/4 to 9/4. Indeed, we have found 
by exact diagonalization that 98 percent of the spectral 
weight of the pinwheel is spread over the energy 0.26 to 
2.50. Also, as shown in Fig. 4, the angle-integrated spec- 
tra at every energy has q dependence that is very close 
to that of isolated dimers. This shows that Delta chains 
have spectra spread over a wide energy range yet have an- 
gle integrated spectral weight essentially that of isolated 
dimers, in agreement with experiments. 

The presence of dilution also provides a simple expla- 
nation for the absence of large peaks at low energy in 
the Raman spectra (37[ that were found in exact diag- 
onalization studies of finite clusters, [l^ The pinning of 
dimers destroys the low energy singlets, by making any 
long-range rearrangement of dimers energetically costly. 
Light scattering should produce two antiparallel triplet 
excitations on neighboring dimers. Two of these spins 
can combine into a singlet to leave a kink-antikink pair 
or two anti-kinks one of which is trapped in a triangle 
with a valence bond. In either case, this leads to energy 
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berg Model is randomly diluted. We have argued that 
this phase has no spin-gap, and supports local two-level 
system excitations as well as singlet and triplet excita- 
tions that decompose into kink-antikink pairs that de- 
localize over a local network of corner sharing triangles. 
This phase provides a consistent picture of thermody- 
namic, neutron and Raman measurements in the Her- 
bertsmithite materials. Significantly reduced concentra- 
tion of impurities can help provide clearer signatures for 
the Valence Bond Crystal phase in the Kagome Lattice 
Heisenberg Model. 

We would like to thank Fabrice Bert, Mark de Vries, 
Joel Helton, Michael Hermele, Andreas Laeuchli, Pe- 
ter Lemmens, Frederic Mila, and Oleg Tchernyshyov for 
many valuable discussions and to Phillipe Mendels for 
an invitation to Paris, where many of these ideas were 
crystalized. 



FIG. 4: Exact diagonalization results for Angle averaged dy- 
namic structure factor for a pinwheel, at different frequencies, 
scaled to have the same maximum, compared with results of 
a single dimer. 



spread over the same range as the triplets, i.e., roughly 
upto 2.5 J. Indeed, the exact diagonalization results of 
Laeuchli and Lhuillier[l3| show that once the low energy 
peaks are removed from the spectra the frequency depen- 
dence of the neutron and Raman scattering profiles are 
very similar. It supports the idea of two weakly interact- 
ing spinous which can form singlets or triplets. Hao and 
Tchernyshyov have shown [15j that these excitations have 
fermionic statistics. The minimum excitation energy is 
reduced when one doesnt have perfect pin-wheel struc- 
tures and the absence of spin-gap means that Raman 
scattering should also extend down to zero energy. These 
features are, indeed, consistent with the experiments. 37 

Since the disordered glass phase breaks no symme- 
try, there should be no phase transition in the system 
and only a gradual crossover from the high temperature 
paramagnet to the VBG phase. The temperature scale 
for the freezing is set by O.IJ, but some strong bonds 
may freeze at temperatures approaching J giving one 
a broad crossover region. Indeed the neutron scatter- 
ing measurements suggest that Dimer correlations per- 
sist beyond T = J/2. On the other hand the Raman 
spectra show a clear development of a quasi-elastic peak 
at low energies as the temperature is raised from 5K to 
295if .[s?! This is suggestive of a gradual melting of the 
Valence Bond Glass into a classical Valence Bond Liquid 
as the frozen dimers arc freed up, giving rise to lots of 
low energy overdamped singlet excitations. Details of the 
Raman spectra, including dependence on polarizationjs^ 
deserve further theoretical attention. 

In conclusion, we have developed a picture for the Va- 
lence Bond Glass phase when the Kagome Lattice Heisen- 



APPENDIX 

In this appendix, we address the recetly reported Neu- 
tron scattering data on the Herbersmithite materials by 
Helton et al. 27| We make the assumptions that the large 
number of free spins created by the impurities go into a 
random singlet phase, where the spins are coupled by a 
renormalized distribution of exchange constants, which 
has a power-law behavior P{J) = CJ'". With this the 
magnetization of the system as a function of field is given 
by 



with y = gh/T with h — ^bH, and 

2 sinh y 



fiiy,x) 



+ 1 + 2 cosh y 
The susceptibility in a field is given by 



with 



X{H,T)T^= / dxx-"f2iy,x) 



f2iy,x) 



2coshy(e^ + 1) + A 
(e^ + 1 + 2coshy)2 



(1) 



(2) 



(3) 



(4) 



Also, the imaginary part of the dynamic susceptibility 
for cj > is given by 



T e"/"^ - 1 



(5) 



Motivated by the experimental observation, we will 
take a — 0.66. We also take g — 2.2. With these the 
plots for M{H, T), x{H, T) and x" are shown in Figures. 
They appear remarkably similar to the experiments. 
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Magnetization versus Field 
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Scaled susceptibility versus h/T 

log-log plot 
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With this assumption, theses spins will contribute 
T^~°' to the specific heat. That combined with a lin- 
ear specific heat coming from the singlet degrees of free- 
dom, can explain the observation of an apparent power- 
law specific heat with an exponent in the range 0.5 to 
0.66.0 



<n/T scaling for dynamical susceptibility 



T 




Experiments show that these power law have a wide 
range of validity in the Herbertsmithite materials. 
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